(3x  mm 

Minm» 


The  University  of  Alberta 
Printing  Department 
Edmonton,  Alberta 


THE  UNIVERSITY  OF  ALBERTA 


THE  PHASE -AMPLITUDE  APPROACH  TO  NONLOCAL  SCATTERING 


by 


WALTER  JOHN  WAT TAMAN I UK 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES 
IN  PARTIAL  FULFILLMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 

OF  MASTER  OF  SCIENCE 


DEPARTMENT  OF  PHYSICS 

EDMONTON,  ALBERTA 


FALL,  1970 


•  • 
11 


UNIVERSITY  OF  ALBERTA 

FACULTY  OF  GRADUATE  STUDIES 

The  undersigned  certify  that  they  have  read,  and 
recommend  to  the  Faculty  of  Graduate  Studies  for 
acceptance,  a  thesis  entitled  THE  PHASE-AMPLITUDE 
APPROACH  TO  NONLOCAL  SCATTERING,  submitted  by  Walter 
John  Wattamaniuk,  in  partial  fulfillment  of  the  require¬ 
ments  for  the  degree  of  Master  of  Science. 


Wvesvi 
lb  P 


ABSTRACT 


The  phase-amplitude  approach  is  considered  as  a 
means  of  calculating  the  regular  and  irregular  wave 
functions  for  nonlocal  scattering.  The  equivalent 
local  wave  function  and  potential  may  then  be  expressed 
exactly  in  terms  of  these.  A  separable  kernel  is  used 
as  an  example.  A  class  of  nonlocal  kernels  that  can  be 
obtained  from  the  Green's  functions  for  second-order 
differential  equations  is  considered.  The  momentum 
dependent  potential  is  derived  and  is  expanded  in  terms 
of  the  index  of  nonlocality  to  obtain  the  effective 
mass  approximation.  A  specific  example,  the  nonlocal 
square  well,  is  solved  analytically  and  its  equivalent 
local  wave  function  and  potential  are  examined  for  the 


s-state . 
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CHAPTER  1 


NONLOCAL  INTERACTIONS 


I .  Introduction 

Historically,  the  use  of  nonlocal  interactions 
began  with  the  advent  of  extensive  research  into  the 
properties  of  nuclear  matter.  Basically,  the  inter¬ 
action  between  complex  nuclei  is  a  many-body  problem 
and  complex  potentials  can  be  calculated,  at  least  in 
principle,  by  summing  over  the  two-body  nucleon-nucleon 
forces.  This  fundamental  approach  has  been  developed 
by  Brueckner ^  and  others and  has  gone  a  long  way 
in  justifying  the  success  of  the  phenomenological  optical 
model  potential  that  is  used  in  the  description  of  nucleon- 
nucleus  interactions.  The  theory  showed  that  additional 
features  should  be  incorporated  into  the  phenomenological 
model.  It  turns  out  that  the  nucleon-nucleus  interaction 
operator  <r|v|r'>,  for  finite  nuclei ^  is  nonlocal  in 
coordinate  space.  This  means  that  the  potential  acting 
on  a  particle  centered  at  a  point  r  depends  not  only  on 
r,  but  also  on  the  value  of  the  wave  function  throughout 
all  the  region  of  nonlocality  and  so  takes  account  of  the 
finite  size  of  the  incident  particle  and  the  dispersive 
nature  of  the  nuclear  matter.  With  the  interaction  now 


. 

' 
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described  by  a  matrix  <r|V|r'>E  V(r,r'),  the  term 
V(r)  ip(r)  in  the  one  particle  Schrodinger  equation 
that  normally  describes  local  interactions  now  becomes 

/  V(r,r 1 )  i|)  (r  ' )  d3r' 

where  V(r,r')  is  the  nonlocal  potential. 

Very  little  is  known  about  the  form  of  this  poten¬ 
tial.  Obviously  it  complicates  matters  because  now  the 
wave  function  satisfies  an  integro-dif f erential  equation 


Brown  and  De  Dominicis 


(3) 


have  shown  that  V(r,r')  must 


be  symmetric  in  r  and  r'  if  the  nucleon  eigen-functions 
are  to  be  orthogonal  to  each  other  within  the  nucleus. 
Furthermore,  the  short  range  of  the  nuclear  forces  and 
the  calculations  of  Brueckner  and  his  coworkers 

/ 

suggest  that  the  range  of  nonlocality  is  small.  If  this 

is  the  case,  then  by  introducing  a  range  of  nonlocality 

as  a  parameter,  one  may  be  able  to  describe  the  effects 
p 

of  the  nonlocality  in  a  phenomenological  treatment  without 
knowing  the  detailed  structure  of  the  interaction  matrix. 
Of  course  one  must  require  that  the  potential  becomes 
local  V (r ,r ' )  =  V(r)  6(r-r')  ,  as  the  nonlocal  range  goes 
to  zero  (3  ->  °°)  . 

Nonlocal  potentials  are  not  confined  only  to  the 
many-body  problem  of  nucleon-nucleus  scattering. 


■ 
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Giltinan  and  Thaler have  demonstrated  that  an 
examination  of  the  energy  dependence  of  the  ^So  and 
"^2  phase  shifts  for  proton-proton  scattering  below 
340  Mev.,  requires  that  the  singlet  even  parity  nucleon- 
nucleon  interaction  be  nonlocal.  This  is  not  surprising, 
for  one  might  expect  that  the  exchange  of  many  particles 
and  resonances  between  two  nucleons  will  produce  non¬ 
local  forces. 

Any  nonlocal  potential  that  is  used  to  describe  a 
scattering  problem  will  have  its  energy  dependent  local 
counterpart.  Indeed,  there  is  no  reason  to  suppose  that 
the  nonlocal  potential  itself  does  not  change  with  energy. 
However  until  such  a  time  as  field  theorists  or  S-matricists 
can  derive  a  nonlocal  potential  from  more  basic  principles, 
a  reasonable  approach  would  be  to  keep  the  nonlocal  poten¬ 
tial  energy  independent  and  vary  it  in  an  appropriate  way 
with  the  least  number  of  parameters  so  that  the  equivalent 
local  potential  changes  with  energy  in  a  natural  way  con¬ 
sistent  with  scattering  data. 

*  *  * 

In  Section  II  we  discuss  the  momentum  equivalence 
of  the  nonlocal  potential  and  obtain  the  nonlocal  radial 
Schrodinger  wave  equation. 

The  phase-amplitude  method  is  considered  in  Chapter 
2  as  a  means  for  obtaining  the  nonlocal  wave  function. 
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We  extend  this  method  in  Chapter  3  to  also  calculate 
the  "irregular"  nonlocal  wave  function.  The  equivalent 
local  wave  function  and  potential  may  then  be  expressed 
exactly  in  terms  of  these.  A  separable  kernel  is  used 
as  an  example. 

A  class  of  nonlocal  kernels  that  can  be  obtained 
from  the  Green's  functions  for  second-order  differential 
equations  is  considered  in  Chapter  4.  The  momentum 
dependent  potential  is  derived  and  is  expanded  in  terms 
of  the  nonlocality  parameter  to  obtain  the  effective  mass 
approximation.  A  specific  example,  the  nonlocal  square 
well,  is  solved  analytically  and  its  equivalent  local 
wave  function  and  potential  are  examined  for  the  s-state. 

II .  Nonlocal  Schrodinger  Wave  Equation 

The  nonlocal  Schrodinger  equation  for  scattering  of 
two  particles  of  equal  mass  M  by  a  nonlocal  potential 
V (r , r ' )  is 

(V2  +  k2)  ip  (r;)  =  f  V  (r ,  r ' )  d3r'  .  (1-1) 

IT 

For  simplicity  we  assume  the  particles  are  spinless  and 
uncharged.  Furthermore  if  we  restrict  ourselves  to 
elastic  scattering,  the  nonlocal  potential  is  real. 

The  wave  function  and  nonlocal  potential  may  be  expanded 


' 


. 
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in  terms  of  spherical  harmonics 


U£(r) 

♦  (e)  =  2  -Jr—Yom(0^) 

Z,m  r  36111 


(1-2) 


M 

ti 


2o+i  K  (r,r')  r.r' 

2  V<E'E'>  =  l  rr'  P£ 


Ko<r'r')  * 

=  rr 1  • 

Jo  f  m 


(1-3) 


The  spherical  harmonics  satisfy  the  differential  equation 


{iIHe  W  sine  h  +  ~TT  r?}  Ytm(0'4>  “ 


sin  0  3(j)' 


-  *U+D  Yfcm(e,<(.) 


and  have  the  orthonormality  conditions 


/  “  ■  «U-  V-  • 

u  u 

Using  these  we  can  substitute  eq.  (1-2)  and  eq.  (1-3)  into 
eq.  (1-1)  to  obtain 

2  oo 

{-^5-  +  k2  -  £(^1)-}u0  (r)  =  /  K.(r,r')  u.  (r')  dr*  . 
dr  r  36  0 

(.1-4 ) 

The  kernel  KA(r,r')  is  calculated  from  V(r,r')  by  inverting 
eq.  (1-3) 


' 
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M 

K^(r,r')  =  2i\rr'  (— ^)  /  V(r,r')  P^(cos0)  d(cos0). 

ft  _1  (1-5) 

Mathematically,  a  nonlocal  potential  is  equivalent 
to  a  momentum  dependent  potential.  This  can  best  be 
illustrated  by  expanding  ^(r1)  about  r'  =  r  in  a 
Taylor's  series. 

ip  (r ' )  =  \p  (r  +  r  '  -  r ) 

=  E  i  [(r 1  -  r).V]ni|;(r) 

XI  • 

n 

=  exp  [(r'  -  r)cV]^(r)  .  (1-6) 

If  we  now  replace  V  by  the  momentum  operator  ^  p  and 
substitute  eq.  (1-6)  into  eq.  (1-1) ,  we  have 

(V2  +  k2)  ip  (r)  =  U(r,p)  ip  (r)  (1-7) 

where  U(r,p)  is  the  momentum  dependent  potential  given 

by 

jj-  (r'-r)  .p  3 

U(r,p)  =  /  V (r ,r 1 )  e  dr  (1-8) 


. 
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CHAPTER  2 

THE  NONLOCAL  PHASE  AND  AMPLITUDE  FUNCTIONS 

I .  Phase  function 

In  Chapter  1,  we  obtained  the  radial  wave  equation 
for  nonlocal  scattering 

2  __  .  00 

{ — 2  +  k  "  ^  }  uo  (r)  =  f  K  (r ,  s)  U  (s)  ds.  (2-1) 

dr  r  36  0  ^  ^ 

We  assume  the  usual  conditions  on  the  nonlocal  kernel 
asymptotically  and  at  the  origin.  These  imply  the 
existence  of  a  regular  solution  u  (r)  characterized  by 

X/ 

the  boundary  conditions 

u£(0)  =  0  •  (2-2) 

U£(r)r  const-sin  [kr  +  ’  (2"3) 

The  constants  6^  are  the  scattering  phase  shifts  and  are 
defined  only  mod  (2tt). 

(5) 

Following  Calogero  ,  we  outline  a  formalism 
whereby  the  phase  shifts  are  calculated  as  the  limits 
of  phase  functions  for  large  r.  For  local  potentials, 
this  method  transforms  the  second  order  differential  wave 
equation  with  two-point  boundary  conditions  into  a  set  of 


. 
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two  coupled  first  order  nonlinear  differential  equations 
with  initial  value  conditions.  This  procedure  greatly 
simplifies  the  study  of  the  analytical  properties  and 
numerical  calculation  of  the  phase  shifts.  For  nonlocal 
potentials  the  generalization  involves  coupled  integral 
equations . 

For  convenience  we  define 


00 

F^(r)  =  /  K£(r,s)  u^(s)  ds 


(2-4) 


With  the  aid  of  the  boundary  condition  at  the  origin  given 
by  eq.  (2-2)  ,  the  nonlocal  wave  equation  may  be  trans¬ 
formed  into  the  integral  equation 


uX.(r)  =  ~  k'  [jj,  (kr)na(ks)  -  nA(kr)  jj,  (ksljF^sJds. 

(2-5) 


The  Riccati-Bessel  functions  j^(z)  and  n^U)  are  defined 
by 


Vz)  =  z  Vz) 

(2-6) 


^(z)  =  Z  n^(z)  . 

where  j  (z)  and  n5  (z)  are  the  spherical  Bessel  functions 
of  the  first  and  second  kind  respectively^  .  We  intro¬ 
duce  the  auxiliary  functions  c£(r)  and  s£(r)  by 


. 
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u£(r)  =  c^(r)  j  ^ (kr )  -  s^  (r)  r\  ^  (kr)  .  (2-7) 

Comparison  with  the  integral  equation  for  u  (r)  implies 

Xj 

1  r  - 

s^(r)  =  -  jr  /  j^(ks)  F^(s)  ds 

(2-8) 

-i  r 

c^(r)  =  1  -  /  n£(ks)  F^(s)  ds 

Thus,  differentiating  with  respect  to  r,  we  have 


si(r)  =  -  I  Vkr>  Vr) 

c'l(r)  =  ■  £ Vkr)  Fjt(r)  ' 

with  the  initial  value  conditions 


(2-9) 


sz(0)  =  0 

(2-10) 

cz(0)  =  1  . 

From  eq.  (2-9)  we  have  the  identity 

c^(r)  j^tkr)  "  s£<r)  ^(kr)  =  0  .  (2-11) 

Using  this  and  the  property  that  the  Riccati-Bessel  func¬ 
tions  satisfy 


(  d2  .  ,2  £  (£+1)  1 

j,(kr) 

(dr*  +  k  -  r2  | 

0£ (kr) 

10 


we  differentiate  eq.  (2-7)  twice  to  find 

u^(r)  +  (k2  -  MA±1L)  u^(r)  =  c^(r)j^(kr)  -  s^  (r )  (kr )  . 

r 

(2-13) 

However,  the  right  hand  side  of  this  expression  is  just 
equal  to  F^(r)  by  comparison  with  eq.  (2-1).  If  we  now 
multiply  eq.  (2-13)  by  u^ (r) ,  we  obtain 


-k[s^  (r)c£  (r)  -  s^(r)c^(r)]  =  F£(r)u£(r)  ,  (2-14) 

where  we  have  made  use  of  eq.  (2-7)  and  the  Wronskian 
identity 


(kr)  (kr)  -  (kr )  (kr )  =  k  . 


(2-15) 


Dividing  eq.  (2-14)  by  u^(r)  yields 


rst(r)i 

i  r 

1 

K(r)J 

k 

S£(r)  - 

Vkr)  -  cTTrT  Vkr) 


F£  (r ) 


Vr) 


(2-16) 


We  define  the  phase  function  6^(r)  by 


sd(r) 

tan  6t( r)  = 


(2-17) 


to  obtain 


Vr)  = 


2  Fit(r> 


-  |  [jJl(kr)cos6Jl(r)  -  (kr) sin6 ^ (r) ]  — ^ 


(2-18) 


. 


, 


- 


subject  to  the  initial  value  condition, 


6^(0)  =  0  . 

The  significance  of  the  phase  function  is  illustrated 
by  expressing  the  wave  function  in  terms  of  it,  that  is 

u*(r)  =  cosl~(rT  l>osVr)  Vkr)  '  sin6Z(r)  Vkr>]  • 

(2-19) 

Since  the  Riccati-Bessel  functions  have  the  asymptotic 
behavior 


5*(kr)  r  oo‘  sin  Ckr  ~ 

(2-20) 

^d(kr)  r  '  cos  ^kr  "  ' 


then  we  have 


c{.(r) 

r  +  »'  cosfi ^  (r) 


sm 


[kr  -  fj-  +  6g  (r)] 


Z 


(2-21) 


If  we  compare  this  with  the  condition  expressed  by  eq. 
(2-3) ,  we  see  that  the  scattering  phase  shifts  are  the 
asymptotic  limits  of  the  phase  function. 


= 


lim 

r+oo 


Vr) 


(2-22) 
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II .  Amplitude  function 

The  amplitude  function  a_  (r)  is  defined  by 

X/ 

a£ (r)  =  (r)/cos6£  (r)  ,  (2-23) 

and  is  more  conveniently  expressed  in  the  form 


a^(r)  =  s^(r)  +  c^(r)  .  (2-24) 

Differentiating  this  equation  and  using  eq.  (2-9)  we 
find 

a^(r)  =  -  — ^ —  [j  £  (kr)  sin6£  (r)  +  n z  (kr)  cos6£  (r)]  . 

(2-25) 

If  we  use  eq.  (2-18) ,  this  can  be  rearranged  in  terms  of 
the  phase  function 


Vr)  =ai 


(r)6^(r) 


j£  (kr) sin6 ^ (r) 
j  ^ (kr) cos6£ (r) 


/S  \ 

+  (kr) cos6^ (r) 

-  (kr) sin6^ (r) j 

(2-26) 


with  the  initial  value  condition 


a£(0)  =  1  .  (2-27) 

The  differential  equations  for  the  phase  and  ampli¬ 
tude  functions  may  be  cast  into  a  more  convenient  form  if 
we  make  use  of  the  Riccati-Bessel  phase  and  amplitude 
functions  D£(z)  and  6^(z)  defined  by 


•1 


1  (?)  ^  -  (1)  tZQOt'dt} 
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3 £  ( 2 )  =  6^  (z)  sinls^  (z) 

(2-28) 

n^U)  =  -  (z)  cosS^(z)  . 

Hence  eqs.  (2-19),  (2-18),  (2-26)  become 


Vr) 

«;<r) 


a£ (r)D^ (kr) sin[6£ (kr) 
“  £  D2 (kr) sin2 [6£ (kr) 


+ 


+ 


Vr)] 

6a(r)] 


(2-29) 

Vr) 

Uj,  <r ) 


(2-30) 


a^(r)  =  -a^(r)5^(r)  cotg[6£(kr)  +  6£(r)]  .  (2-31) 


For  the  s-state  (£=0) ,  the  above  equations  simplify  to 
u  (r)  =  a(r)  sin  [kr  +  6 (r)  ] 

6 '  (r)  =  -  i  sin2  [kr  +  6  (r)  ]  (2-32) 

a'  (r)  =  -  a(r)  cotg  [kr  +  6(r)]  . 

Furthermore,  if  the  potential  is  local 

i.e.  K(r,s)  =  V(r)  6 (r-s)  f 

then  the  s-state  phase  function  satisfies  the  first-order 
nonlinear  differential  equation 

6'(r)  =  -  V^-—  sin2  [kr  +  6(r)]  (2-33) 


with  6(0)  =  0 . 


' 
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Analytical  solutions  to  the  above  equation  (and  in 
general  for  any  £  wave)  are  of  course  confined  to  poten¬ 
tials  for  which  the  Schrodinger  equation  is  a  priori 
solvable.  However , numerical  solutions  may  now  be  obtained 
in  a  direct  and  rapid  manner  for  any  physical  potential. 


III.  Integro-dif f erential  equation  for  6^(r) 

If  the  nonlocal  phase  function  is  known,  the  ampli¬ 
tude  function  is  determined  from  eq.  (2-31)  by  integration. 
Using  the  condition  a^(0)=l,  we  have 

r 

a^(r)  =  exp  {-  /  6  £ (s) cotg  [6 ^ (ks)  +  6^(s)]ds}  .  (2-34) 

Using  the  definition  of  F^(r)  in  eq.  (2-4)  we  now  have 


pyr) 

“*(r) 


/  ds  K^(r,s)D^(ks) sin 


.exp  {  /  6^ (t) cotg  [6 0 (kt) 


£ 


(ks)  +  6 ^  ( s) ]  . 
+  6£(t)]dt}  . 


(2-35) 


The  integrand  of  the  exponential  term  is  not  singular, 
because  the  poles  of  cotg[6jl(kt)  +  S£(t)]  are  cancelled 
by  zeros  of  6j(t).  The  phase  function  may  now  be 
expressed  as  the  solution  of  the  integro-dif f erential 
equation 


. 


' 
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6I(r)  =  "  k  sin[6£(kr)  +  6  £  (r )  ] 

oo 

•  /  K£(r,s)  D£(ks)  sin[6£(ks)  +  S£(s)]ds 

r 

•  exp  {  /  6£(t)  cotg[6£(kt)  +  6  (t)]dt}  .  (2-36) 


IV.  Born  and  improved  Born  Approximation 

Except  for  special  forms  of  the  nonlocal  kernel, 
the  integro-dif f erential  equation  is  probably  impossible 
to  solve.  One  must  resort  to  approximation  methods  to 
evaluate  the  phase  function  for  arbitrary  potential.  The 
most  straightforward  of  these  is  the  method  of  successive 
approximations.  That  is;  the  successively  better  appro¬ 
ximations  are  given  by  6»n+^  (r)  where 


6^(n+1)(r)  =  -  |S£(kr)  sin[6£(kr) 

OO 

./  K£  (r ,  s)  D£  (ks)  sin  [<5 £  (ks ) 

. exp{  f  6£n^  (t)  cotg [6 £ (kt) 
s 

To  start  the  iteration  we  use  6  ^ 
we  have 


+  «Jn) (r)] 

+  (s) ] ds 

+  S^n) (t) ] dt} .  (2-37) 

(r)  =  0.  To  first  order 


fij(r)  ~  _  I  ^(kr)  Kj,(r’s)  ds 


(2-38) 


. 


. 
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and  the  phase  shifts  are  therefore 

0°  °0  ^ 

<$n  ~  “  it  /  /  K  (r ,  s)  j  (kr)  j  (ks)  dr  ds  .  (2-39) 

0  0  36  36 

However,  these  are  just  the  phase  shifts  in  the  Born 
approximation,  for  if  we  consider  the  eq.  (2-16)  for 
the  tangent  of  the  phase  function 

i  2 

[tamS^fr)]  =  -  jj-  [j^tkr)  -  tan6£ (r)  (kr) ] 


Vr> 

U£(r) 

(2-40) 


and  take  as  zero-order  approximation  6 „  ( r )  =  0,  then  to 

X/ 

first  order 


tan6„  = 


k  0 


00 

/ 

0 


K 


£ 


(r ,  s)  j  ^  (kr )  j  ^  (ks)  dr  ds 


(2-41) 


This  reassures  us  that  the  validity  of  the  Born  approxi¬ 
mation  requires  that  the  phase  shifts  be  small,  so  that 
the  tangent  does  not  differ  appreciably  from  its  argument. 

In  any  case,  one  should  be  able  to  show  that  the 
iteration  scheme  outlined  for  the  phase  function  converges 
to  the  correct  solution  for  any  physical  kernel.  However, 
because  of  the  strong  nonlinearity  of  the  phase  equation, 
the  resulting  expressions  already  become  unmanageably 
complicated  by  the  second  cycle.  One  must  ultimately 


' 
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resort  to  numerical  calculation  to  use  the  successive 
approximations  scheme. 

The  Born  approximation  can  be  improved  upon  by 
making  the  approximations 


[j^  (kr) -tan6^(r)  (kr )  ]  2  ~  j^(kr)-2  tan6^  (r)  (kr)  j  Q  (kr) 


(2-42) 


U£ (r )  ~  j^(kr)  , 


thereby  effectively  linearizing  the  equation  for  tan60(r) 


[tan6£(r)]  ~  [:p(kr)  -  2  tan6 p (r ) n p (kr) ] 


k  LJ£ 


£ 


00 


./  K£(r,s)  j  ^  (ks)  ds 


(2-43) 


The  solution  to  this  equation  is  given  by 


00 


tan6^(r)  = 


1  r  /s 

-  V  f  j o (ks) ds/K  (s,t) j p (kt)  dt 

K  o  *  0 


00 


.e 


p  r  - 

£  /  n^(kx)dx  /  K£(x,y) j^(ky)  dy 
s  0 


(2-44) 


and  the  resulting  phase  shifts  are  determined  by  taking 


the  limit. 


■ 

, 
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V.  Separable  kernel 

One  example  of  a  nonlocal  potential  where  the  phase 
and  amplitude  functions  can  be  determined  analytically  is 
that  of  a  separable  kernel, 

K£(r,s)  =  g£(r)  g£(s) 

From  eq.  (2-4)  we  have 

Vr)  =  g*(r) 

where  the  constants  are  defined  by, 

00 

=  k //  g£(s)  u£(s)  ds 

Eqs .  (2-8),  (2-17)  and  (2-24)  then  yield 

r  „ 

-  /  j£(ks)  g£(s)  ds 

tan6^(r)  =  — - — - - 

A£  -  /  n£(ks)  g&(s)  ds 

and 

«!.<*>  = 1  -  ^  l  Vks)  g*(s)  ds  + 

r  o  r  ~  2 

(/  n0 (ks)g0 (s)ds)  +  (/  j  (ks)g.(s)ds) 

0  ^  0 


(2-45) 


(2-46) 


(2-47) 


(2-48) 


(2-49) 


' 


. 
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The  constants  A^  are  obtained  by  making  use  of  the 
integral  equation  for  the  wave  function  (eq.  (2-5)) in 
eq.  (2-41) 

°o  s  ^ 

k+ f  g  (s)  ds  /  {  j  (ks)  r|  (kt) -n  (ks)  j  „  (kt)  }g  (t)  dt 

qX/  q  X/  ^  X/  X/  X/ 

a  =  - 

Z  oo 

/  g^(s)  j^(ks)  ds 


(2-50) 
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CHAPTER  3 


EQUIVALENT  LOCAL  POTENTIALS 


The  use  of  nonlocal  potentials  in  nuclear  physics 
makes  a  detailed  study  of  their  connection  with  local 
potentials  worthwhile.  We  define  an  equivalent  local 
potential  (ELP)  by  the  requirement  that  it  reproduces 
the  same  asymptotic  wave  function  as  the  nonlocal  poten¬ 
tial.  Mathematically  we  require  real  functions  x^(r) 
and  (r)  that  satisfy  the  equations 


{— +  k2  -  (r)  =  V  (r)  x^r) 

dr  r 


(3-1) 


x*(r)  F 


CO 


sin 


[kr  -  . 


(3-2) 


Here  Vn (r)  is  the  ELP  to  be  determined,  x  (r)  is  the 
regular  equivalent  local  wave  function  and  6^  are  the 
phase  shifts  determined  from  the  asymptotic  behavior  of 
the  nonlocal  wave  function  u^(r). 

It  is  immediately  obvious  from  eq.  (2-1)  that 
trivial  solutions  are 


x*(r) 

Vr> 


const . u^ (r) 


f  K. (r , s) u  (s)  ds 
0  * _ _ _ 

u*(r) 


(3-3) 


. 


. 
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However,  in  this  case,  (r)  has  a  pole  whenever  u^(r) 

has  a  zero  and  is  therefore  far  from  satisfactory. 

Ideally  we  require  an  ELP  that  has  no  poles  (unless 

perhaps  at  the  origin)  and  falls  away  rapidly  for  a 

short  range  kernel.  We  now  outline  an  analysis  due  to 
( 7 ) 

Fiedeldey  for  obtaining  ELP ' s  that  satisfy  these 

requirements,  although  they  will  be  implicit  functions 
of  the  scattering  energy  and  angular  momentum. 


I .  Equivalent  Local  Potentials  expressed  in  terms  of 
nonlocal  wave  functions 

We  assume  that  there  exist  two  linearly  independent 
solutions  of  the  equivalent  local  wave  equation  (3-1) ,  a 
regular  solution  x^(r)  satisfying  eq.  (3-2),  and  an 
irregular  solution  y^(r),  that  has  the  asymptotic  behavior 


y„  (r)  - *  cos  [kr  -  4^-  +  60]  (3-4) 

2  £,  r  00  L  2  jcJ 

and  whose  "local  Wronskian"  with  x^(r)  is 

x^(r)  y£(r)  -  x^(r)y^(r)  =  k  .  (3-5) 

We  shall  relate  the  nonlocal  wave  function  u^ (r)  to  the 
equivalent  local  one  by  an  interpolating  function  f £ (r) , 

u£(r)  =  f£(r)  x£(r)  . 


(3-6) 


' 


* 
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Substituting  eq.  (3-6)  into  the  nonlocal  wave  equation 
(2-1) ,  we  find 


{-%  +  k2  -  Ml+il}  x  (r)  = 

dr  r ^  1 


00 


f— Yr)  (r ,s) (s)x^ (s)ds-f^ (r)x£ (r) -2f ^ (r)x^ (r) } , 

0 


(3-7) 


and  from  eq.  (3-1) ,  the  ELP  must  be 


oo 


Vr)  = 


f  K  (r , s) f  (s) x  (s)ds-f " (r)x  (r) -2f ! (r)x' (r) 

q  X/  X/  X/  X/  Xs  X/X/ 


f*(r)  Vr) 


(3-8) 


With  the  aid  of  the  local  Wronskian  (eq.  (3-5) ) ,  we  intro¬ 
duce  the  identity 


kx„(s)  =  [x’(r)  y„(s)  -  x£(s)  yj(r)]x„(r) 


l '  'J  l 


+  [yj,(r)  x^fs)  -  x^tr)  y^(s)]  x^(r) 


(3-9) 


The  ELP  then  becomes 


Vr)  =  FTrT  {.  Vr'S)  Vs)  [X£  (r)  y^  (s)  -Xj  (s)  y^  (r)  Jds-f  ^  (r) 

X  *  (  IT  )  oo 

+  FTrT  (ST(7r){^  f0  K4(r,s)fJl(s)[yJl(r)xJl(s)-xJl(r)yJl(s)]ds-2f;( 


(3-10) 


' 


■ 


. 
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We  now  try  to  choose  f^(r)  in  such  a  way  that  there  are 
no  singularities  in  the  ELP .  If  we  assume  that  f  (r)>0, 
then  the  absence  of  poles  is  ensured  by  requiring  that 


fs.(r)  = 


_1_ 

2k 


/  (r,  s)  f  g  (s)  [y0  (r)xQ  (s)  -  x9  (r)  y  Q  (s)  ]  ds  , 


0 


Z 


z 


z 


z 


(3-11) 


and  it  follows  that 
f - (r)  , 

VJl(r)  =  ~  ft(r)  +  kf4(r) 

oo 

./  K£(r,s)f  ^(s)  [x^(r)y£(s)-x£(s)y^(r)]ds  .  (3-12) 

In  the  limiting  case  of  a  local  kernel,  (r , s) =V (r ) 6 (r-s) , 

eqs ,  (3-11)  and  (3-12)  must  reduce  to 

f  (r)  =  constant 

Xj 

(3-13) 

V£(r)  =  V (r) 

With  the  aid  of  eq.(3-5)  it  is  easily  seen  that  these  condi¬ 
tions  are  fulfilled.  Equation  (3-11)  is  particularly 
useful  in  establishing  that  the  function  (r)  defined  by 

Vr)=Vr>  y*(r) 


(3-14) 


•  • 
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also  satisfies  the  nonlocal  wave  equation.  This  may  be 
easily  proven  by  considering  the  "nonlocal  Wronskian" 

u^(r)  v£(r)  -  u£(r)  v^(r)  =  k  f^(r)  (3-15) 

which  is  derived  from  the  definitions  supplied  by  eqs. 
(3-6)  and  (3-14)  and  with  the  aid  of  eq.  (3-5) .  If  we 
differentiate  eq.  (3-15)  and  use  eq.  (3-11)  to  replace 

i 

f^(r)  it  follows  that 
u^(r)  v£(r)  -  u£(r)  v£(r)  = 

oo  OO 

/  K^(rfs)u^(s)ds  -  u^(r)  /  (r,s)v^(s)ds  .  (3-16) 

If  we  now  use  the  nonlocal  wave  equation  (eq.  (2-1))  to 
replace  u^(r),  it  follows  that  v^(r)  satisfies  the  equa- 
tion 

00 

v”(r)  +  {k2  -  =  /  K{, (r  's)  v  sl  (s)  ds  - 

(3-17) 

However,  because  of  the  normalization  imposed  on  y^(r) 
by  eq.  (3-5) ,  the  "irregular"  solution  to  the  nonlocal 
wave  equation,  v^(r) ,  will  have  the  asymptotic  behavior 

v^fr)  V00)  cos  Ekr  "  X  +  5J' 


(3-18) 


' 


. 
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where  and  6^  are  now  determined  from  the  asymptotic 

behavior  of  the  regular  nonlocal  wave  function  u^(r),  i.e. 


u«.(r)  r  ->•  »"  f«,(”)  sin  [kr 


(3-19) 


Assuming  that  we  have  a  priori  determined  u^(r) 
and  hence  v^(r)  according  to  the  prescription  above,  the 
interpolating  function,  f^(r)  and  the  ELP,  V^(r),  can 
now  be  obtained  exactly  in  terms  of  them.  f^(r)  i-s  °f 
course  determined  from  the  nonlocal  Wronskian  (eq.  (3-15)). 
The  expression  for  the  ELP  may  be  simplified  using  the 
equation 


,  u*tr)  Vr)  Vr) 

x£(r)  =  f  (r)' - 2 - 

*  Vr)  f£(r) 


and  a  similar  equation  for  y^(r).  We  have  then 


Vr)  - 


f^(r) 

f£(r) 


oo 


kf £  ^ 


/  K£(r,s)  [u^(r)v^(s)-v^(r)u£(s)]ds 


0 


l£ 


£ 


£ 


kf2(r) 


( c'  t  \ 

^f£(r)l 

f,(D 


/  K  (r,s)  [v^trJUj,  (s) -u£  (r)  v£  (s)]  ds  . 

(3-20) 


However, from  eq.  (3—16),  we  note  that  the  integral  in  the 
above  equation  is  just  the  differential  of  the  nonlocal 
Wronskian  (eq.  (3-15)).  Thus  we  have 


. 


i  &[_(  -  '<*)  '/-(f  )  <  )  a  ($'>*!  ‘ 


. 
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Vr)  = 


f;<r) 

f  )l  <r) 


+  2 


Vr) 


00 


kf2(r)  0  1 


f  KQ  (r,  s)  [u'Q  (r)  (s)  -vj  (r)  uQ  (s)]  ds  . 


Z 


(3-21) 


A  simpler  expression  for  the  ELP  can  be  obtained 
by  considering,  as  solutions  for  the  equivalent  local 
wave  equation,  the  following, 


x^(r)  =  A^(r)  sin  B^(r) 
y^(r)  =  A^(r)  cos  B^(r) 


where 


As. (r) 


■£  00 


1 


Vr) 


- *•  kr 

r  ■>  oo 


(3-22) 


(3-23) 


If  we  now  substitute  for  x^(r)  i-n  ec3*  (3-1)  we  obtain 
for  the  ELP , 


V4(r)  =  k2  -  ’ 


AI(r) 

vn 


-  [Bl(r)]2  + 


+  Vr) 


2A^  (r) 

Vr) 


Bj(r)l 


cotg  B^(r) 


(3-24) 


A  similar  ELP  is  obtained  for  y£(r)  except  that  cotg  B£(r) 
is  replaced  by  -  tan  B  (r)  .  It  is  obvious  that  to  have 
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both  x^(r)  and  y^(r)  satisfying  the  same  equation,  and 
to  avoid  untold  numbers  of  singularities  due  to  the 
tangent  and  cotangent  functions,  we  must  have 


2  A^(r)  Bj(r) 

Vr>  +  = 


(3-25) 


However, this  condition  is  already  implied  by  the  Wronskian 
prescription  in  eq.  (3-5).  That  is,  we  have 

A^(r)  B^(r)  =  k  (3-26) 

and  eq.  (3-25)  follows  by  differentiation.  The  ELP 
becomes 


Vr) 


Z  ( £+1 ) 


AI(r) 

Vr) 


Al(r) 


(3-27) 


where,  from  eq.  (3-22) 


A'(r) 


xl{  r)  +  yj(r) 


In  terms  of  the  nonlocal  wave  functions  we  have 


A^(r) 


U£(r)  +  VliT) 


(3-28) 


It  is  easily  verified  that  with  A£(r)  as  defined  above, 
the  ELP  given  by  eq.  (3-21)  reduces  to  eq.  (3-27)  as  it 
must.  Equation  (3-27)  for  the  ELP  is  not  very  useful 
physically  since  it  seems  to  indicate  that  the  potential 


■ 


. 


. 

. 
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depends  explicitly  on  the  scattering  energy  and  angular 
momentum.  This  is  not  the  case,  for  eqs.  (3-15)  and 
(3-21)  for  f^(r)  and  V^(r),  show  that  the  range  for  which 


they  differ  from  a  constant  and  zero  respectively,  coin¬ 
cides  with  that  of  the  nonlocal  potential. 

II .  The  irregular  nonlocal  wave  function 

Since  we  require  the  irregular  nonlocal  wave  function 
to  calculate  the  ELP,  we  shall  now  show  how  Calogero's 
phase-amplitude  method  can  be  extended  to  calculate  v^(r). 


The  integral  equation  for  v^(r)  can  be  written  as 


.G^(s)  ds 


(3-29) 


where  a0  and  b  are  constants  to  be  determined  and  G„(r) 

Xj  X/  ^ 

is  defined  by 


CO 


G£(r)  =  f  K£(r,s)  v^(s)  ds  . 


(3-30) 


As  in  Chapter  2,  we  define  the  auxiliary  functions 


a 


b^  -  £  /  j£(ks)  ds 


(3-31) 


' 


V 
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so  that  eq.  (3-29)  may  now  be  written  as 


v£(r)  =  -  c0(r)n0(kr)  -  s0(r)j0(kr) 


l 


(3-32) 


The  functions  s£(r)  and  c£(r)  here  defined  have  no  connec¬ 
tion  with  the  ones  used  in  Chapter  2.  Since  their  use  is 
only  auxiliary,  we  shall  have  no  conflict  of  definitions. 
From  eqs.  (3-31)  we  have 


s^r)  =  |  n£(kr)  G£(r) 


c£(r)  = 


- 1  Vkr>  Vr) 


(3-33) 


and  hence  the  identity 


s£  (r )  j£(kr)  +  c^r)  ^(kr)  =  0  ’  (3-34) 

Using  this  identity  and  eq.  (2-12) ,  we  differentiate  eq. 
(3-32)  twice  to  find 


v"(r)  +  {k2  -  us(r)  =  ”  ct(r)^«.(kr)  “  si(r)^S,(kr)  ' 

36  r 

(3-35) 

and  therefore  by  comparison  with  eq.  (3-17) 

G£ (r)  =  -  c£(r)  n^(kr)  -  s£(r)  j£(kr)  .  (3-36) 

With  the  aid  of  eqs.  (2-15),  (3-32),  and  (3-34),  the  above 

equation  may  be  written  as 


■ 


' 
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s  (r)  • 

k(W’  = 


sji(r)  - 

;Vkr)  -  <?7(7r)  Vkr) 


n2  G  (r) 


Vr) 


.  (3-37) 


Defining  the  "irregular  phase  function"  y  (r)  by 

X/ 


s  (r) 

tan  Vr)  =  E7TFT  ' 


(3-38) 


eq.  (3-37)  becomes 


Y;<r)  = 


1  r-  -  -,2  Go<r) 

-  ^  [njl(kr)cosY£(r)  +  jA(kr)sinYjl(r)]  — m 

^  Xj 


(3-39) 


The  irregular  wave  function  is  now  written  as 


Vr)  = 


Vr) 


cosYt(r)  [-n^(kr)cosY)l(r)  -  j£  (kr)  sinY^  (r)]  , 

(3-40) 


and  hence  we  define  the  "irregular  amplitude  function" 
S£(r)  by 

eA(r)  =  c^(r)/cosy£(r)  , 
which  is  more  conveniently  written  as 

$£ (r )  =  c^(r)  +  s^(r)  .  (3-41) 

We  differentiate  this  expression  and  use  eq.  (3-33)  to 
obtain 


BJ,(r)  =  £  tn^(kr)sin6j,  (r)  -  j ^  (kr)  cos6^  (r)]  (r)  .  (3-42) 


■ 
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Eq.  (3-39)  is  now  used  to  replace  G£(r)  in  eq.  (3-42) 


=3£(r) Y^(r) 


j£  (kr)cosy£ (r)  -  f|£  (kr)  siny£  (r ) 

^  /N 

,-n^ (kr) cosy£ (r)  -  j ^ (kr) siny  (r) 


(3-43) 

We  shall  now  make  use  of  the  Riccati-Bessel  phase  and 
amplitude  functions  defined  by  eq.  (2-28)  to  rewrite 
eqs.  (3-39),  (3-40),  and  (3-43). 


v£(r)  =  $£(r)D£(kr)  cos[<S£(kr)  +  y£(r)] 

G  ( r ) 

=  -  !S'<kr>  ^2[6£(kr)  +  Yjl(r)](^J-7r) 

X/ 

3^(r)  =  e^(r)y^(r)  tan[6£(kr)  +  y£(r)]  . 


(3-44) 


(3-45) 


(3-46) 


These  equations  are  very  similar  in  form  to  their  counter¬ 
parts  (eqs.  (2-29)  to  (2-31))  obtained  for  the  regular 
nonlocal  wave  function.  There  is  the  important  difference 
that  instead  of  initial  value  conditions  as  for  a^(r)  and 
6^(r)  we  now  have  end  value  conditions  on  $^(r)  and  y^(r). 
To  see  this  more  clearly  we  make  use  of  the  asymptotic 

/N  A 

expansions  of  6^(kr)  and  (kr)  which  may  be  determined 
from  eqs.  (2-20)  and  (2-28). 


' 


D£(kr) 


r  ->  oo 


1 


Vkr)  kr 


£tt 

2 


(3-47) 


Therefore  the  irregular  wave  function  has  the  asymptotic 
behavior 

V!l(r)  r  +  =o‘  S£(“)  cos  [kr  -  ^  +  y£(<»)]  (3-48) 

and  according  to  eqs.  (3—18)  and  (3—19)  we  must  therefore 
have  the  end  value  conditions 


6,(r) 


- > 

r  oo 


a£  (°°) 


Yj,(r) 


r  oo 


(3-49) 


where 


a£  (°°) 


lim  a 

2T-3-OQ 


£ 


(r) 


6 


£ 


lim  6£(r) 


(3-50) 


If  one  attempts  a  numerical  calculation,  the  sequence  of 
steps  is  quite  clear.  One  integrates  the  equations 
( (2 . 29) - (2 . 31) )  to  some  point, r  say,  where  the  kernel 
is  effectively  zero  and  obtains  a£(r)  and  S£(r).  Then 
using  the  end  value  conditions  3£(r)  =  a£(r),  and 
y  (r)  =  60(r),  one  integrates  the  equations  ( (3-44) - (3 . 46) ) 
backwards  to  the  origin  to  evaluate  vp (r) . 


■* 
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Because  of  the  structural  similarity  of  the  phase- 
amplitude  equations  for  the  regular  and  irregular  non¬ 
local  wave  functions  one  can  obtain  a  similar  form  of 


integro-dif f erential  equation  for  y  (r)  and  proceed  to 

X/ 


discuss  Born  and  improved  Born  approximations  for  tany  (r) . 

Xj 


However,  the  generalizations  are  quite  straightforward  and 
no  attempt  shall  be  made  here  to  deal  with  them. 

Ill .  The  nonlocal  Wronskian  expressed  in  terms  of  phase- 
amplitude  functions 

If  we  differentiate  the  phase  amplitude  expression 
for  u£(r)  given  by  eq.  (2-29) ,  we  have 


u£(r)  =  (a£(r)D£(kr)  +  ct£  (r)  D£  (kr )  }  sin  [6  £  (kr)  +  6£(r)] 

+  a£(r)D£(kr) (6£(kr)  +  6£ (r) }cos  [S£ (kr)  +  6£(r)]. 

(3-51) 


(2-31)  to  express  a' (r)  in  terms  of  6'(r),  we 


Using  eq. 


find  that  the  terms  containing  S£(r)  cancel  and  we  can 


write 


D£(kr) 


u£(r)+a£(r)D£(kr)  6£  (kr)  cos  [6£  (kr)  +  6£  (r)]  . 

(3-52) 


In  a  similar  fashion  we  differentiate  eq.  (3-44)  for  v£(r) 
and  use  eq.  (3-46)  to  obtain 


■ 


- 
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•  BI(kr) 

Vr)  =  fT^T)  vt(r)_Pt(r,Dt(kr,8i(kr)sin[«Jl(kr)+YJl(r)]. 

(3-53) 

The  nonlocal  Wronskian  then  becomes 

U*(r)  Vr)-U£(rW(r)=Vr)  Vr^^krJ^tkrJcos^  (r)-Y£  (r)]  . 

(3-54) 

From  eq .  (2—28)  it  may  easily  be  verified  that 

^£(kr)  ^(kr)  =  k  '  (3-55) 

therefore 

ui (r)  vi (r)  ~us. (r)  vi (r)  =ka£ (r)  ai (r)  cos  [6s, tr)  ~y<i (r)]  • 

(3-56) 

It  follows  from  eq.  (3-15)  that  the  interpolating  function 
f^(r)  is  given  by  the  simpler  expression, 

f  ^  (r )  =  a^(r)3£(r)  cos[<S£(r)  -  y£(r)]  .  (3-57) 

Since  the  nature  of  the  differential  equations  that  a  (r) 

A/ 

and  3 ^ ( r )  satisfy  ensure  that  both  are  positive  and 
greater  than  zero,  we  conclude  that  Fiedeldey's  method 
of  obtaining  an  ELP  is  valid  only  for  nonlocal  potentials 
such  that 


1 <$£(r)  ”  Y£(r)  I  <  7  • 


' 


» 


' 
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If  this  condition  is  violated  then  we  may  expect  poles 
to  appear  in  the  ELP .  We  must  also  ensure  that  f  (r) 

X/ 

has  the  correct  asymptotic  behavior.  If  we  take  the 
limit  of  large  r  and  use  eg.  (3-49) ,we  find 


lim  f  (r ) 

r+oo 


2  ,  » 


which  is  the  correct  normalization  required  by  eg.  (3-19) 


IV.  Interpolating  function  and  ELP  for  a  separable  kernel 

For  the  separable  kernel  given  by  eq.  (2-45) ,  we  can 
immediately  write  down  the  solution  to  the  irregular  wave 
function  from  eg.  (3-29) . 


VJ>(r)  =  a^SL(kr)  ~  bJtrU(kr)  " 


-  jp  /[]A(kr)Tif  (ksM^ksln^krjJg^sJds 


£ 


£ 


£ 


£ 


(3-60) 


where  the  constants  B.  are  defined  by 


oo 


B 


/  9^(r)v£(r)  dr 


(3-61) 


According  to  the  prescription  given  by  eqs.  (3  31) ,  (3—38) , 

and  (3-41)  we  have 


-®*V  f0  Vks>Vs)as 


b£B£  -  f  j£(ks)g£(s)ds 


tany^ (r) 


(3-62) 


- 


and 


»'( r>= 


1  r  - 

V  f0  3  s,  (ks)  (s)  ds 


2  2T  2 

+  f'V  ^  £  n^ksjg^sjds1 


(3-63) 


For  convenience  we  define  the  functions 


s0 (r)  =  /  j  o (ks)g£ (s)  ds  , 


0 


(3-64) 


ci  (r>  = 


-  /  n£(ks)  g£(s)  ds 


We  now  impose  the  end  value  conditions  in  eg.  (3-49) . 
Thus , referring  to  eqs.  (2-48)  and  (3-62)  and  imposing 
the  condition,  tan6^(°°)  =  tany^M,  we  have 


ajBt  +  Cl("} 

V”’  -  biBl 


s 


+ 


l 


(o°) 

cl(”) 


(3-65) 


The  constant  B „  is  determined  using  eqs.  (3-60)  and  (3-61) 


BAS*(~>  +  b£Ct(”>]  =  k  +  fQ  g£(S)  dS 

.  /  [j  ^  (ks)  (kt)  -  j  ^  (kt)  ri^  (ks)  ]  (t)  dt 

However  using  eq.  (2-50)  this  may  be  written  as 


The  nonlocal  Wronskian  (eq.  (3-15))  may  be  calculated 
directly  by  using  eq.  (2-5)  and  eq.  (3-60)  for  the  non¬ 
local  wave  functions.  The  result  is 


2  SJ.(r)  afSf.(r)  +  b8C8.(r) 

_  T_  _  Ay  |  Ay  A/  As  Ay 

JT  '  "  £  bT“  An 


(3-67) 


and  thus 


f>)  = 

Making  use  of  eq. 


aHS8,(”)  +  bAC  <CO> 

— +  rt  • 

(3-66)  in  eq.  (3-67)  we  have  finally 


f£(<»)  =  bs.  •  (3-68) 

2  2 

Now  we  make  use  of  the  end  value  condition,  f  (°°)  =a ^  (°°)  , 
to  obtain 


2,  > 

2  C£  °° 
b  =  1  +  —  c  (oo)  +  - 2 

£  A£  £  a 


+  s^(°°) 


£ 


(3-69) 


The  above  equation  follows  from  eq.  (2-49)  for  a^(r) . 
Substituting  for  b^  into  eqs.  (3-66)  and  (3-67)  we  also 

have 


c  m 

'iVr1 

S£V  1 


s?(oo)  +  cP(oo) 

( _£ _ ± _ ) 

s  (») 


(3-70) 


and 


. 
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i 

a*=A*  V“> 


)  (c^  (°°)  +S^  (°°)  ) 


SA(00) 


(3-71) 


The  interpolating  function  and  the  ELP  are  then  given 
by 


9  s p  (°°)  c  (r)  -  c  (°°)  s  (r) 

f*<r>  =  b*  (1  +  -  A^(~) - — > 


and 


f0(r>  f»(r)  2  b0  r  -  ,  1 

Vfc(r)=-f^TrT  +  + 

gz(r) 

•  ■!<-> 
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CHAPTER  4 


ANALYSIS  OF  A  NONLOCAL  KERNEL 


We  consider  the  class  of  nonlocal  kernels  that 
satisfy  the  differential  equation 

{  TT  "  "  ~U21)  }  KP  (r,r')w"i5(r)  =  -X3V5  (r  ' )  6  (r-r  ' 


(4-1) 


with  homogeneous  boundary  conditions 


K£(0,r')  =  K£(«,  r')  =  0 


(4-2) 


The  solutions  to  eq.  (4-1)  are  given  by 

K£  (r ' r ' )  =  A3wi5(r)w^(r,)3£(i3r<)h^1)  (i3r>)  ,  (4-3) 

where  w(r)  is  a  real  short  range  function  characteristic 
of  a  physical  kernel  and  h^1^  (z)  is  the  Riccati-Hankel 
function  of  the  first  kind  defined  by 

(z)  =  Vz)  +  1  ^ (z)  • 

This  example  is  a  special  case  of  a  more  general 

( 8 ) 

class  of  nonlocal  kernels  considered  by  Razavyv  '  that 
can  be  obtained  from  the  Green  functions  of  second  order 
differential  equations.  The  kernel  is  physically 
realistic  in  the  sense  that  the  parameter  X  measures 


' 


1 
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the  strength  of  the  interaction  while  —  is  a  measure 

p 

of  the  range  of  nonlocality  of  the  potential.  That  is; 
if  we  take  the  limit  as  jjj-  0,  we  have  from  eq.  (4-1) 


K£(r,r')  =  Aw  (r ) 5 (r-r ' ) 


and  the  interaction  becomes  local. 

The  choice  of  the  "shape  function"  w(r)  is  res¬ 
tricted  by  the  boundary  conditions  imposed  upon  the 
kernel  (eq.  (4-2)).  In  addition,  we  must  ensure  the 
physical  convergence  of  the  phase-amplitude  equations, 
(eqs.  (2-30),  (2-31)).  Using  the  expansions  of  the 

Riccati-Bessel  functions  asymptotically  and  at  the  origin, 
it  is  easy  to  show  that  these  requirements  are  met  if  we 

demand  that  w(r)  diverge  less  than  — at  the  origin  and 

rz 

vanish  faster  than  —  at  infinity. 


I .  Nonlocal  and  momentum  dependent  potential 

The  nonlocal  potential  V(r,r')  is  easily  obtained 

by  substituting  the  expression  for  the  kernel  (eq.(4-3)) 

(9) 

into  eq.  (1-3)  and  using  the  expansion 


3  j  r-r ' |  3  rr '  &=0 


00  *  ^  n  )  r .  r 

- —  I  (2)1+1)  j t(iBr<)h^i)  (ierjp,  (-^rr)  . 


>  r  rr 


We  have 


M 


—o  V(r  ,r ' ) 
fT 


w*5  (r)  w*3  (r ' )  £ 
4tt 


-3|r-r ' 


(4-4) 


r-r 


■ 
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Since  the  Yukawa  term  in  eq.  (4-4)  satisfies 


=  -4tt6  (r-r '  ) 


r-r ' 


then 


h 


or  formally 


The  nonlocal  Schrodinger  equation  (eq.  (1-1))  becomes 


(V2  +  k2)  i|i  (r) 


and  the  momentum  dependent  potential  defined  in  eq.  (1-7) 
is 

2  i 

t  P  “1  v 

U  (r  ,p)  =  A  w 2  (r)  (1  +  ■  9---9)  w2(r) 

3  fi 

II .  The  effective  mass  approximation 

If  one  assumes  that  the  nonlocality  is  small,  that 
2  2  2 

is,  3  >>|Aw(r) |  and  3  >>  k  ,  then  one  may  expand  the 
operator  (1  -  V2/32)  1  to  first  order  in  1/32 


2  -1 


V 


2 


1  + 


2 


3 


V 


■ 
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Therefore  eq.  (4-5)  becomes 
2  2 

(V  +  k  )ij>(r)  =  X  W(r)ij>(r)  +  (r)  V2w^  (r )  (r )  .  (4-6) 

e 

With  the  aid  of  the  operator  identity  (here  primes  refer 
to  differentiation  with  respect  to  r) 


k  2  k  k  ?  ..  [w'  (r)l 2 

w^(r){V  w2(r)-w2(r)V  }  =  ^{w  (r) - }  + 


2w  (r ) 


+  w'  (r)  (i  +  JL) 
r  9r 


eq.  (4-6)  may  be  rewritten  as 


V.  {  (1  -  — y  — )  Vip  (r)  }  +  k2i|i  (r) 

6 

l  [w '  ( r )  ]  2 
=  A{w(r)  "  —2  L  2w(r) 

Z  P 


11  /  \ 

-  w  (r)  - 


2w ' (r) 


■]  }ip(r)  . 


(4-7) 


The  above  equation  may  be  regarded  as  a  generalization 


of  the  ordinary  local  Schrodinger  equation  with  the  mass 
M  replaced  by  an  effective  mass 


M*  = 


M 


1  - 


X  w(r) 

e2 


which  is  a  function  of  position.  Eq.  (4-7)  becomes 
-*2(7  M*7F7  7)*(E)  +  ^{w(r)  ‘  ^2^W2WrT-  "  W"(r)- 


_  2w'  (r)j  j  ^  (r )  =  E  i(i(r) 


(4-8) 


, 


' 
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This  is  the  well  known  effective  mass  approximation  of 
Frahn  and  Lemmer  . 


Ill .  Equivalent  local  potential  in  the  effective  mass 
approximation 


It  is  possible  to  obtain  an  ELP  in  the  effective 
mass  approximation  by  making  the  transformation 

i(i  (r)  =  (1  -  Xw^r)  )  *  <p(r)  .  (4-9) 

8 

Here,  <p  (r)  is  the  equivalent  local  wave  function,  for, 
if  we  use  eq.  (4-9)  directly  in  eq.  (4-7) ,  we  find  that 
i|;(r)  satisfies  the  Schrodinger  equation 


V2cj)  (r)  +  {k2-  V (r)  }  <j)  (r)  =  0  . 


V(r)  is  the  energy  dependent  ELP  given  by 


V  (r)  =  X{  w  (r)  (1  -  K)  (1  -  ^L^L) 

3  3^ 


[W  (r)] 
4^2w(r) 


-}(1  - 


Xw  (r ) x  2 

2  } 

3 

(4-10) 


The  difference  Xw(r)— V(r)  may  be  expressed  as 

.  ,  2  [«’ (r)] 2 

Xw(r)  -V  (r)  =  (1  -  ^^(k2-  Xw(r))  +  2  ,T}  x 

qz  3^  43  w (r ) 

Xw (r ) * “2 

2 

3 


x  (1 


(4-11) 


. 


. 
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We  note  that  the  ELP  becomes  w(r)  when  the  range  of  non¬ 
locality  (g)  is  zero  (as  required)  .  Since  w(r)  is  a 
positive  shape  function  the  right  hand  side  of  eq.  (4-11) 
is  negative  for  attractive  potentials  (X  <  0) .  In  this 
case,  the  behavior  of  eq.  (4-11)  implies  that  for  fixed 
w(r) ,  the  average  strength  of  the  interaction  decreases 
as  a)  the  range  of  nonlocality  (g)  increases  and  b)  the 
scattering  energy  increases. 


IV.  The  nonlocal  square  well 

If  one  multiplies  the  radial  equation  for  nonlocal 
scattering  (eq.  (2-1))  by  w  2(r),  the  differential  equa¬ 
tion  for  the  kernel  (eq.  (4-1))  may  be  used  to  obtain  a 
fourth  order  differential  equation  for  the  nonlocal  wave 
function 

{-4  -  e2  -  +  k2  -  Uji(r) 

dr  r  dr  r 

=  -  X32w^ (r)u^ (r)  .  (4-12 

One  of  the  few  cases  for  which  analytic  solutions 
to  eq.  (4-13)  exist  is  that  of  the  square  well  of  range 

b 

1  r  <  b 


w  (r) 


0 


r  >  b 


' 


. 
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Since  the  local  potential  well  is  known  to  be  a  rough 
approximation  to  the  actual  interaction  between  nucleons, 
it  is  worthwhile  to  study  this  more  sophisticated  non¬ 
local  model.  De  La  Cruz  et  al^11^  attempted  to  fit  the 
experimental  phase  shifts  for  proton-proton  scattering 
in  the  dS  state  using  an  attractive  nonlocal  square  well 
in  a  range  intermediate  between  a  repulsive  hard  core  and 
an  attractive  Eckart  potential.  They  deduced  a  set  of 
parameters  that  was  effective  in  fitting  the  phase  shifts 
for  laboratory  energies  in  the  range  20-340  Mev. 

We  shall  study  only  the  pure  nonlocal  square  well 
here  with  the  specific  purpose  of  determining  its  equi¬ 
valent  local  wave  function  and  potential.  For  the  s- 
state,  eq.  (4-12)  reduces  to 

2  2 

(--d -  -  £2)  +  k2)  u  (r )  =  -X32u(r)  ,  (r  <  b)  .  (4-13) 

dr  dr 

The  solutions  to  eq.  (4-13)  are  expressed  as  linear  com¬ 
binations  of  the  functions  sinhyr,  coshyr,  sinvr  and 
cosvr,  where 

u2  =  >5  (e2-  k2  +  J(g2+k2) 2  -  4g2x  } 


v2  =  h  {k2 


e2  +  J(e2+  k2)2  -  4g2x  }. 


(4-14) 


.  pH " 
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For  attractive  potentials  (A.  <  0)  ,  y^  and  are  always 
positive . 

The  wave  function  is  uniquely  determined  by  the 
boundary  conditions 


u(r)  F^Tf 

kr 

u" (0)  +  k2 

u(0) 

=  0 

u”(b)  +  k2 

u(b) 

=  0 

where  the  last  two  equations  are  a  consequence  of  eq. (4-2) . 
The  wave  function  is  written  as 


u(r)  = 


k  sinvb  sinhyb 


2  2  2  2 
v(y  +k  ) sinhyb  +  y (v  -  k  ) sinvb 


x 


x 


sinhyr 

sinhyb 


2  2 
(U  +  JO 


sinvr -i 
sinvo 


(4-15) 


The  phase  and  amplitude  functions  are  then  immediately 
determined  from  the  equations 

ot(r)  sinfkr  +  S(r)]  =  u(r) 

k  ctn  [kr  +  6(r)]  =  d^r~  • 

The  irregular  solution  to  the  nonlocal  wave  func¬ 
tion,  which  we  require  in  order  to  compute  the  equivalent 
quantities,  is  given  by 


!>4*  *  >  %;?53Klr 


»• 

■  1  ftJb 


■ 
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v(r)  -  A  sinhyr  +  B  coshyr  +  C  sinvr  +  D  cosvr  . 


(4-16) 

The  boundary  conditions 

v"  (0)  +  k2  v(0)  =  0 
v"  (b)  +  k2  v(b)  =  0 

imply 


A 


,cosvb-coshybN „ 

C"  sinhyb  )B  + 


( 


2 

v  - 


y2+ 


sinvb 

sinhyb 


(4-17) 


The  constants  B  and  C  are  determined  according  to  the 
prescription  given  in  Section  II  of  Chapter  3.  That  is, 
we  determine  the  irregular  phase  and  amplitude  functions 
from  the  equations 


B(r)  cos  [kr  +  y(r)]  =  v(r) 

-k  tan  [kr  +  y(r)]  = 

The  end  value  conditions,  3(t>)=a(b),  and  y(b)=6(b),  then 
determine  the  constants  B  and  C.  However,  it  is  unneces¬ 
sary  to  do  this  calculation,  for  if  we  directly  use 
eqs.  (4-15) ,  (4-16)  and  (4-17)  to  evaluate  the  nonlocal 


Wronskian,  we  obtain 


. 
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2 

kf  (r)  =  u ' (r ) v (r )  -  u(r)v'(r)  (4-18) 

2  2 

=  kB(^-+v2)  {1  +  — _ _ (v2~k2)  (y2+k2)T(r) _ } 

v  “k  (H  )  [v  (y2+k2)  sinhyb+y  (v2-k2)  sinvb] 

where  the  function  T(r)  is  given  by 


T (r )  =  v (sinhyr  cosv (b-r ) +cosvr  sinhy (b-r ) -sinhyb) 

(4-19) 

+  y (sinvr  coshy (b-r ) +coshyr  sinv (b-r ) -sinvb) 


Since  T(0)  =  T(b)  =0,  we  have 

f2(o>  =  f2(b>  =  B 

v  -kz 

and  from  eq.  (4-14) 

(v2  -  k2)  (p2  +  k2)  =  -xe2  . 
Therefore  eq.  (4-18)  can  be  written  as 


(4-20) 


,f (r) , 2_ 
^ (b) '  ~ 


1  - 


A3  T (r) 


(y  +v  )[v(y  +k  ) sinhyb+y (v-k  ) sinvb] 


(4-21) 


V.  Properties  of  the  interpolating  function  f (r) 

The  function  T(r)  defined  by  eq.  (4-19)  is  symmetri¬ 
cal  about  r=b/2.  Also  the  upper  and  lower  bounds  on  T(r) 


' 


# 
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are  quite  easily  obtained.  We  can  rewrite  eq.  (4-19)  in 
two  different  ways 

T(r)  =  (vsinhyr-ysinvr) (cosv (b-r) -coshy (b-r) ) 


+  (vsinhy (b— r) — ysinv (b-r) ) (cosvr-coshy r ) 


(4-22) 


and 


T (r) +2 (vsinhyb+y sinvb) = (vsinhyr+y sinvr ) (cosv (b-r) +coshy (b-r) ) 


+  (vsinhy (b-r ) +ysinv (b-r )) (cosvr+coshyr )  .  (4-23) 


Equations  (4-22)  and  (4-23)  follow  from  eq.  (4-19)  upon 
replacing  b  by  r+b-r  and  expanding  sinvb  and  sinhyb. 
Because  of  the  inequalities 


vsinhyx  :>,  ysinvx 
coshyx  ^  cosvx 


(0  <  x  ^  b)  , 


we  can  see  by  inspection  that  the  right  hand  side  of 
eq.  (4-22)  is  always  negative,  while  the  right  hand  side 
of  eq.  (4-23)  is  always  positive.  It  follows  that 

-2 (vsinhyb  +  ysinvb)  ^  T(r)  ^  0  .  (4-24) 

Eq.  (4-24)  can  be  used  directly  in  eq.  (4-21)  to  obtain 
bounds  on  the  ratio  f(r)/f(b).  With  the  aid  of  eq.(4-14) 


we  have 


H 


' 
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2  2  2  2.  29 

/B  +k  ^  v (y  +k  ) sinhyb-y (v  -k  ) sinvb 

v  2  2  *  2  2 - 2 2 - 

y  +v  v(y  +k  ) sinhyb+y (v-k) sinvb 


(jLrJ..)  2 

'f  (b)  j 


(4-25) 

It  is  convenient  to  define  a  renormalized  nonlocal  wave 
function  u(r)  by  u (r ) =u (r) /f (b) .  Since  u (r) =f (r) x (r) 
we  have 


u(r)  =  f  (r)  _ _ 

x  (r)  f  (b)  r  -►  b 


(4-26) 


There  are  two  observations  of  importance  that  can 

be  made  regarding  eq.  (4-25) .  The  first  is  the  fact  that 

2 

the  lower  limit  of  (f(r)/f(b))  is  always  greater  than 
zero  and  therefore  no  poles  can  appear  in  the  ELP .  For 
our  particular  example,  Fiedeldey's  method  of  constructing 
an  ELP  is  entirely  consistent  and  the  assumption  that  f (r) 
be  greater  than  zero  is  not  necessary. 

The  second  observation  is  that  the  ratio  f(r)/f(b)= 
u(r)/x(r)  is  always  less  than  unity  within  the  potential 
well  and  therefore  the  nonlocal  wave  function  u(r)  is 
always  smaller  than  the  equivalent  local  wave  function 
x(r)  there.  This  result  supports  an  argument  in  a  paper 
by  Austern^^  where  an  approximate  procedure  is  used  to 
construct  the  equivalent  local  wave  function.  He  points 
out  that,  in  general,  an  attractive  nonlocal  potential 
equivalent  to  a  local  potential  (in  the  sense  that  they 


' 


, 
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give  rise  to  the  same  asymptotic  wave  functions  for  all 
energies)  produces  a  wave  function  that  is  always  smaller 
inside  the  region  of  the  potential  than  the  equivalent 
local  wave  function. 

VI .  The  ELP  for  the  nonlocal  square  well 

The  equivalent  local  potential  for  the  square  well 
is  determined  from  the  expression 

V (r)  =  -  £1^  +  2(0^-)2  +  — i—  /  K(r,s)  [u-  (r)v(s)  - 
f  (r)  kr  (r )  0 

-  v  '  (r)  u  (s)  ]  ds 


Because  u(r)  and  v(r)  both  satisfy  the  nonlocal  wave 
equation  (eq.  (2-1)),  the  integral  in  the  above  equation 
can  be  rewritten  in  terms  of  differentials  of  u(r)  and 
v (r) .  Making  use  of  the  nonlocal  Wronskian  (eq.  (3-15)) , 

we  obtain 


V  (r) 


f"  (r) 
f  (r) 


+  2  ( 


f'  (rK 

f  (rP 


2 


u1  (r)  v"  (r)  -v  1  (r)  u"  (r)  +  k2 
kf2(r) 


(4-27) 


The  first  two  terms  in  the  above  equation  are  evaluated 

(4-21) .  The  last  two  terms  are 


by  differentiating  eq. 


' 


' 
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evaluated  using  eqs.  (4-15) ,  (4-16) ,  (4-17)  and  (4-20) . 
After  a  somewhat  tedious  calculation  we  obtain 


V  (r) 


=  3 


2  (yS2 (r) -vS1 (r) )  3  4  S^(r) 


2Q(r) 


+  J  ^ 


Q2  (r) 


(4-28) 


— 2 - ^ - {v  (k2+y 2 )  (sinhyb-S-,  (r)  )  +y  (k2-v2)  (sinvb-S-  (r)  }  , 

(y  +v  )Q(r) 


where  for  convenience  we  have  defined 


Q  (r )  =  (|-^y)  2  [v  (y2+k2)  sinhyb+y  (v2-k2)  sinvb]  , 

S1(r)  =  sinhyr  cosv(b-r)  +  cosvr  sinhy (b-r)  , 

(4-29) 

S  (r)  =  sinvr  coshy (b-r)  +  coshyr  sinv(b-r)  , 

S3(r)  =  sinhyr  sinv(b-r)  -  sinvr  sinhy (b-r)  . 

2 

For  large  6  ( 62>>k2 , | X | ) ,  the  expressions  for  p 

2 

and  v2  in  eq.  (4-14)  may  be  expanded  in  powers  of  1/3  • 
To  lowest  order  we  have 

2  _  a2  i 

y  £  3  -  A 

v2  ~  k2  -  A  , 

and  from  eq.  (4-22)  we  see  that 


T  (r)  x  -/k2  -  X  sinh3b  . 


. 
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Therefore,  the  ratio  of  nonlocal  to  local  wave  functions 
becomes 


(u(r)_ 
lx(r)  ; 


(b)  ; 


~  1  + 


3- 


However ,  this  ratio  is  equal  to  that  obtained  in  the 
effective  mass  approximation  (eq.  (4-9)) ,  where  we  had 


( 


ijj  (r)  2 


4>  (r) 


■)  = 


(1  -  -\) 
3 


=  1  +  At 
3 


To  evaluate  the  ELP  for  large  3/  we  note  that  the 

first  term  in  eq.  (4-28)  is  proportional  to  T" (r)  while 

the  second  is  proportional  to  JjT’ (r)]  .  Because  the 

function  T(r)  becomes  a  constant  at  large  3,  these  terms 

effectively  vanish.  Specifically,  it  can  be  shown  that 

2 

they  vanish  faster  than  1/3  for  large  3.  For  the 
remaining  term  in  eq.  (4-28)  we  can  write 

vtk^+y^) (sinhyb-S^ (r) +p (k^-v2) (sinvb-S2 (r) ) 


=  -k  T(r)-yv[(ysinhyr+vsinvr) (cosv (b-r) -coshy (b-r) ) 

+  (vsinv (b-r) tysinhy (b-r) ) (cosvr-coshyr)] 
zjk2-  X  (32+  k2-  X)  sinh3b  . 


Now,  since 


.  I 


' 


■ 


2 

Q  (r )  s  3 2  ( 1  +  -4-)  (1  +  h-ZL) 


3 


3' 


sinh3b 


t 


we  have 


V  (r )  ~ 


(1  +  —■)  (1  +  —J21) 


»  X(1  + 

3 


3‘ 


3' 


This  is  just  the  equivalent  local  potential  in  the  effec¬ 
tive  mass  approximation  because  from  eq.  (4-10)  we  have 


[V  (r )  ] 


ef f .mass . approx . 


X(1  -  fx> 


(  - 


-L) 

3 


~  X  (1  + 


l-k‘ 

32 


Therefore ,  within  the  range  of  the  interaction  (0  <  r  <  b) 
we  have  verified  that  the  ELP  obtained  from  Fiedeldey's 
method  matches  that  of  the  effective  mass  approximation. 

The  ELP  for  the  nonlocal  square  well  (eq.  (4-28)) 
is  a  rather  formidable  function  of  radial  distance, 
scattering  energy,  nonlocality  and  well  depth  and  range. 

No  attempt  is  made  here  to  study  its  behavior  analytically 
A  numerical  calculation  of  the  ratio  V(r)/\  was  made  for 
the  purpose  of  extracting  some  of  the  qualitative  features 
of  the  ELP.  The  range  of  values  studied  was 

1  <  E  <256  Mev. 

-640  <  v  <  -10  MeV. 
o 

%  <  $  <  16  (fermi)  1 

1  <  b  <  5  fermi  , 


■ 

' 

I ' , 
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where  vq=  Xfi  /M,  E  =  and  M  is  the  nucleon  mass. 

Figures  (1.-3.)  offer  in  part  some  idea  of  the  behavior 
of  the  ELP.  It  was  generally  found  that  the  ratio 
V(r)/X 

a)  is  always  less  than  unity  and  decreases  as  the 
range  of  nonlocality  (1/3)  increases  (Fig.  1.), 

b)  decreases  as  the  scattering  energy  increases 

(Fig.  2. )  , 

c)  decreases  as  the  well  depth  | X |  increases 
(Fig .  3  . )  . 

An  obvious  disadvantage  of  the  equivalent  local 
potentials  in  figures  (1.-3.)  is  the  fact  that  they  have 
cusps  of  finite  height  and  width  at  r  ~  0  and  r  b. 

This  irregular  behavior  is  a  consequence .  of  our  choosing 
as  a  shape  function  w(r),  the. unit  step  function,  which 
has  a  discontinuity  at  r=b.  This  discontinuity  makes  the 
ratio  f"(r)/f(r)  in  eq.  (4-27)  to  be  finite  at  r=0  and 
r=b  and  thus  causes  the  cusp.  Of  course,  as  the  inter¬ 
action  becomes  local  (1/3  0)  t  the  width  of  the  cusp 

vanishes  and  the  step  function  is  retained  as  the  ELP . 
This  suggests  that  Fiedeldey's  method  should  more  appre¬ 
ciably  be  applied  to  shape  functions  that  decay  smoothly 
to  zero.  Unfortunately  it  would  be  extremely  difficult 
to  analyse  such  nonlocal  potentials  analytically . 


o 
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